Soliton turbulences in the complex Ginzburg-Landau equation 



Hidetsugu Sakaguchi 
Department of Applied Science for Electronics and Materials, 
Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 
Kasuga, Fukuoka 816-8580, Japan 

We study spatio-temporal chaos in the complex Ginzburg-Landau equation in parameter regions 
of weak amplification and viscosity. Turbulent states involving many soliton-like pulses appear 
in the parameter range, because the complex Ginzburg-Landau equation is close to the nonlinear 
Schrodinger equation. We find that the distributions of amplitude and wavenumber of pulses depend 
only on the ratio of the two parameters of the amplification and the viscosity. This implies that 
a one-parameter family of soliton turbulence states characterized by different distributions of the 
soliton parameters exists continuously around the completely integrable system. 
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Dissipative structures and chaos have been studied in nonlinear- nonequilibrium systems. Various kinds of spatio- 
temporal chaos or weak turbulences were found in Rayleigh-Benard convections, electro-hydrodynamic convections of 
liquid crystals and chemical systems [H, 0. The Kuramoto-Sivashinsky equation is one of the simplest model which 
exhibits spatio-temporal chaos, and the statistical properties of the spatio-temporal chaos have been intensively 
studied [|[ 0, [|J M . The complex Ginzburg-Landau equation has been also intensively studied as a model of spatio- 
temporal chaos [3, [H, 0|. The Kuramoto-Sivashinsky equation appears as a phase equation for the complex Ginzburg- 
p ^ Landau equation near the onset of the phase instability. The energy spectrum of the spatio-temporal chaos in 
the Kuramoto-Sivashinsky equation exhibits almost flat in a large scale (small wavenumber regime), which has an 
analogy with the equi-partition law of energy in thermal equilibrium states. On the other hand, in hydrodynamics 
. turbulences, which are described by the Navier-Stokes equation, the energy spectrum has a singular form called the 

1 " i ' Kolmogorov spectrum and the strong turbulences exhibit multi-fractal intermittency • These behaviors are closely 

related to the singular behavior of the Euler equation, which appears in the limit of no viscosity of the Navier- 
Stokes equation. In the limit of no viscosity and no amplification, the complex Ginzburg-Landau equation is reduced 
to the nonlinear Schrodinger equation, which is well known as a completely integral system and a typical soliton 
equation. In this brief report, we study the statistical properties of the complex Ginzburg-Landau equation in the 
regime of weak amplification and viscosity. The complex Ginzburg-Landau equation with weak amplification and 
viscosity appears naturally in some problems of nonlinear optics such as optical fibers with dissipation and external 
pumping 11]. Kishiba et al. tried to explain the energy spectrum in the complex Ginzburg-Landau equation with 
weak amplification, saturation and viscosity [l2j • In the spatio-temporal chaos in the regime of weak amplification and 
viscosity, soliton-like pulses play an important role and it is rather different from the spatio-temporal chaos studied 
before in Refs.[8] and [9]. We will show its peculiar statistical properties of the soliton turbulences and discuss a 
L" ■ relation with the behavior in the nonlinear Schrodinger equation. 
. £^ ' Our model equation has the form 
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where is a complex variable, ei and C2 are respectively parameters of amplification and viscosity. We consider the 
case of £i > and £2 > in this paper. In the limit of ei = £2 = 0, Eq. (1) is reduced to the nonlinear Schrodinger 
equation and has a family of soliton solutions: 

A Ak(x — vt)—iujt 
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where A is the amplitude of the soliton, k is the wavenumber, and v = k is the velocity of the soliton. The amplitude 
and the wavenumber can be independently changed as soliton parameters. The nonlinear Schrodinger equation is a 
completely integrable system for infinite length (L = 00). Then, the nonlinear Schrodinger equation has an infinite 
number of invariants including the total norm ./V = /_ |</>| 2 da; and the total momentum P = J_ 00 (id4>* /dxcf> — 
id4>/dx<p*)/2dx, and so on. If e\ > or e2 > 0, Eq. (1) becomes a dissipative system, and the invariants of motion 
disappear. 

We have performed numerical simulations of Eq. (1) with the split-step Fourier method. The timestep is At = 0.001, 
and periodic boundary conditions are imposed. The system size L is a control parameter. Even for very small 
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FIG. 1: (a) Time evolution of \<j>\ at e\ — £2 = 0.01 and L = 6. (b) Time evolution of peak amplitudes |0(x)| at £1 = £2 = 0.01 
and L = 6. The solid (dashed) line denotes the peak amplitude of the pulse located in x < L/2 (x > L/2). (c) Time evolution 
of peak positions of the two pulses for L — 6.7 at ei = £2 = 0.01. 
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FIG. 2: (a) Time evolution of |<jf>| for L = 20 at £1 = £2 = 0.01. (b) Snapshot profile of \<f>(x, t)\ at t = 2150. Dashed curves 
denote approximation by soliton-type solutions, (c) Time evolution of \<f>\ for L = 20 at £1 = £2 = 0.. 



parameters of t\ and £2, the chaotic behavior appears for large system. Regular behavior appears in a small system. 
Figure 1(a) displays a temporal evolution of |0| at ej = £2 = 0.01 and L = 6. Two pulses appear from a uniform state 
= 0, because = is an unstable solution. The peak positions do not change in time at L = 6. However, the peak 
amplitude of the pulses is breathing in time. The breathing motions of the two pulses are out of phase as in seen in 
Figure 1(b). As L is increased, the peak positions begin to move. An example of moving pulses is shown in Fig. 1(c) 
at L = 6.7. The spatial motion of the two pulses is almost synchronized but the breathing motions are out of phase 
in time. 

When the system size L is further increased, more pulses are created and spatio-temporal chaos appears. Figure 2(a) 
displays chaotic time evolution of \4>{x, t)\ for L = 20 at £1 = £2 = 0.01. Creation and annihilation of pulse structures 
occur, and radiation-like waves with small amplitudes also appear. Figure 2(b) is a snapshot profile of \<f>(x, t)\ at t = 
2150. Three dashed curves denote 2.3/ cosh{2.3(x- 6.32)}, 2.95/ cosh{2.95(x- 11.9)} and 1.76/ cosh{1.76(x- 17.47)}. 
The pulse structures are well approximated by the soliton solutions Eq. (2), because the parameters £1 and £2 are 
rather small. The spatiotemporal chaos might be therefore interpreted as soliton turbulence. Figure 2(c) displays 
time evolution of \<p(x, t)\ for t > 2000 in another numerical simulation, in which t\ = £2 takes the same value 0.01 
as in the case of Fig. 2(a) before t = 1900, but £1 = £2 = are set to zero after t = 1900. That is, <j> obeys the 
nonlinear Schrodinger equation after t = 1900. The initial condition at the fresh start time t = 1900 is a state which 
has appeared as a result of the spatio-temporal chaos. The time evolution in Fig. 2(c) is not chaotic, because the 
nonlinear Schrodinger equation is a completely integrable system. However, the time evolution of \4>\ in Fig. 2(c) looks 
similar to Fig. 2(a). 

To characterize statistical properties of the spatio-temporal chaos, we have calculated the modulus A = \(f>\ at the 
local maxmum points in the profile \<f>(x, t)\, which is interpreted as an amplitude of soliton, and a local wavenumber 
of soliton: k — (— i/2)(0 x 0* — 4>* x 4>) / \<t>\ 2 at t = 5n (n is an integer). And we have constructed a distribution P(A) 
and P(k) of A and k in a larger system L — 80. The parameters £1 and £ 2 are changed as £1 — £2 = 0.01, 0.001, 0.0001 
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FIG. 3: (a) Distributions P(A) of amplitude of pulses, (b) Distributions P(k) of wavenumber of pulses. Time evolution of 
|0| for L = 20 at ei = £2 = 0.01. Eight distributions for eight parameter sets of £i = £2 = 0.01,0.001,0.0001 and 0, and 
ei = £2/2 = 0.01,0.001,0.0001 and are plotted. The ratio r = £i/£2 for the parameter sets is 1 or 1/2, and the value is 
denoted near the distributions. 
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FIG. 4: (a) Average value of A as a function of £1 for £2 = £i and £2 = 2e\. (b) Root mean square of wavenumber k as a 
function of ei for £2 = £1 and £2 = 2ei . (c) The first Lyapunov exponent as a function of £1 for £2 = £1 and £2 = 2ei . 



and 0, and t\ — £2/2 = 0.01,0.001,0.0001 and 0. The two initial conditions in the case of t\ = £2 = are the 
two final states which were numerically obtained for the parameters £1 = £2 = 0.0001 and e\ = £2/2 = 0.0001. 
The numerical simulation was performed until t = 30000 and the data between t — 20000 and 30000 were used to 
construct the distributions. Figures 3(a) and (b) display the distribution P(A) and P(k) for the eight parameter 
sets. The distributions P(A) and P(k) overlap very well, when the ratio r = ei/e 2 takes the same value 1 or 
1/2. The overlap is seen even in the limit of t\ = £2 = 0. This is consistent with the result seen in Figs. 2(a) 
and (c). That is, the statistical properties of the spatio-temporal chaos and the integrable system are almost the 
same, if the initial condition for the nonlinear Schrodinger equation is chosen as a state in the spatio-temporal 
chaos. On the other hand, the distributions for ei/e2 = 1 and 1/2 are definitely different. We have calculated the 
average value of A and the root mean square of k for e x = e 2 = 0.01,0.005,0.0025,0.001,0.0005,0.00025,0.0001,0 
and £1 = £ 2 /2 = 0.01,0.005,0.0025,0.001,0.0005,0.00025,0.0001,0. The results are shown in Figs. 4(a) and (b). 
Almost horizontal lines imply that the distributions of P(A) and P(k) depend only on the ratio r = £i/f2, when £1 
is sufficiently small. We have also calculated the first Lyapunov exponent which characterizes the spatio-temporal 
chaos. Small deviation 5<fi from <fi(x,t) obeys a linearized equation of Eq. (1): 

.w) iw) / a 2 wn r _. 

The first Lyapunov exponent was calculated as the average value of the linear growth rate of the quantity 
[Jo Id^dx] 1 / 2 . Figure 4(c) displays the first Lyapunov exponent for the same parameter sets as in Figs. 4(a) and 
(b). The first Lyapunov exponent increases from as a function of ei, which implies that the spatio-temporal chaos 
becomes stronger as £1 is increased. The first Lyapunov exponent is naturally for t\ = £2 = 0. These behaviors seem 
to be strange, but they are not paradoxical, because the irregularity is not generated from a regular initial condition 
owing to the non-positive Lyapunov exponent, but the regularity is neither created from an irregular initial condition 
owing to the non-negative Lyapunov exponent in the nonlinear Schrodinger equation. 

The distributions P(A) and P(k) depend on the ratio r — ei/ea continuously. Figures 5(a) and 5(b) display the 
average value of A and the root mean square of k as a function of the ratio r = ei/e2 for a fixed value of £2 = 0.001 . The 
average value of A and the root mean square of /c is a increasing function of the ratio r. The total norm N = \4>\ 2 dx 
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FIG. 5: (a) Average value of A as a function of the ration r = ei/e2 for 62 = 0.001. The dashed curve is (3r) 1/2 . (b) Root 
mean square of wavenumber k as a function of the ratio r = 61/62 for 62 = 0.001. 
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— = [{d(j)/dt)(t>* + 4>{d(t>/dt)*]dx = / {2e 1 \cj)\ 2 -2e 2 \d(j)/dx\ 2 )dx. (4) 
Jo Jo 

If <p(x, t) is approximated by a sum of solitons as 4>(x) ~ Y^i Ai exp{iki(x— Xi)}sech{A(x— Xi)}, Eq. (4) is approximated 
as 

^- = 2 J2(eiA i -e 2 Al/3-e 2 k 2 A i ). (5) 

i 

If the total norm is assumed to be constant in time, the distances \x% — Xj\ between two solitons are large, (^4 3 ) 
is approximated by {A) 3 and (fc 2 ) is neglected, (A) is estimated as {A) ~ (3ei/e 2 ) 1 ^ 2 = (37-) 1 / 2 , which is a rough 
estimate of the pulse amplitude and denoted by the dashed curve in Fig. 5(a). We have not succeeded in estimating 
the root mean square of k, yet. 

To summarize, we have found that soliton turbulence appears when ei and e 2 are sufficiently small in the complex 
Ginzburg-Landau equation (1). The soliton turbulence is a chaotic attractor of the complex Ginzburg-Landau equa- 
tion, so the time evolution leads to the chaotic attractor from almost all initial conditions. The soliton turbulence is 
characterized by definite distributions of P(A) and P(k), and the distributions are determined by the ratio r = e\/e 2 - 
Even if t\ is decreased to zero with a fixed ratio r — e\/e 2 , the distributions P(A) and P(k) take almost the same 
form. On the other hand, the Lyapunov exponent is decreased to zero, when ei is deacreased to zero. In the limit 
of ei = and e 2 = 0, the complex Ginzburg-Landau equation is reduced to the nonlinear Schrodinger equation. In 
the nonlinear Schrodinger equation, there is no attractor and no ergodicity, and the time evolution is completely 
determined by the initial conditions. That is, a one-parameter family of soliton turbulence states corresponding to 
the different ratio r = £i/e 2 , which are characterized by different distributions of A and k, exists around the com- 
pletely integrable system. This is a unique relation between the weak spatio-temporal chaos exhibited by the complex 
Ginzburg-Landau equation and the completely integrable dynamics exhibited by the nonlinear Schrodinger equation. 
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